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Abstract 

We investigate tunneling properties of a bound pair of Fermi atoms in an optical lattice, com- 
paring with results obtained in an attractive Hubbard model. In the strong coupling regime of 
the Hubbard model, it has been predicted that the motion of a bound pair between lattice sites 
is accompanied by virtual dissociation. To explore the possibility of this interesting phenomenon 
in optical lattice, we calculate molecular wavefunction in a cosine-shape periodic potential. We 
show that the molecular tunneling accompanied by dissociation occurs in the intermediate coupling 
regime of the optical lattice system. In the strong coupling regime, in contrast to the prediction 
in the Hubbard model, the bound pair is shown to tunnel through lattice potential without disso- 
ciation. As a result, the magnitude of molecular band mass M remains finite even in the strong 
coupling limit, which is in contrast to the diverging molecular mass in the case of the Hubbard 
model. Including this finite value of molecular band mass, we evaluate the superfluid phase tran- 
sition temperature T c in the BEC limit of the optical lattice system, where the Hubbard model 
gives T c = due to the diverging molecular mass. 

PACS numbers: 03.75.Ss, 71.10.Ca, 37.10.Jk 



I. INTRODUCTION 



In 2006, the superfluid state was realized in a 6 Li Fermi gas loaded on an optical 



ii. 



lattice [1|, l2j. In the optical lattice, atoms feel a periodic potential produced by standing 
wave of laser light [3 L Thus, together with a tunable pairing interaction associated with a 
Feshbach resonance , we can now study lattice effects on Fermi superffuids in the 

BCS-BEC crossover region. Since Fermi gases in optical lattices are similar to conduction 
electrons in metals, superfluid lattice Fermi gases may be also useful for the study of metallic 
sup er conduct ivity. 



Various effects of optica 



theoretically [7, 
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attice on superfluid Fermi gases have been studied 
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14 



15] . Among them, one important effect is an 



anisotropic Fermi surface. When the lattice potential is strong, the lattice Fermi gas is ex- 
pected to be close to the Hubbard model, consisting of nearest-neighbor hopping —t and on- 
site pairing interaction —U. In this case, the Fermi surface in the cubic lattice has the nesting 
property at the half-filling, characterized by the nesting vector Q = (7r/d,7r/d,7r/d), where 
d is the lattice constant. This perfect nesting induces strong density wave fluctuations, the 
strength of which is comparable to pairing fluctuations [l2j, [l6j. As a result, the coexistence 
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16] . 



of superfluid state and density wave state is realized in the half-filling casejll. 
Since the competition between these two kinds of fluctuations is absent in a uniform Fermi 
superfluid, this coexistence phenomenon is characteristic of lattice Fermi superffuids. We 
note that the anisotropic Fermi surface has been recently observed in a 40 K lattice Fermi 



gas 



13]. 



Besides this, kinetic properties of bound pairs (Cooper pairs) are also strongly affected 
by optical lattice potential. In the strong coupling regime of the Hubbard model, it has 
been shown that the hopping of a bound pair between lattice sites is accompanied by virtual 
dissociation This comes from the fact that the ordinary Hubbard model consists of 
the atomic hopping and on-site interaction. Namely, when the bound pair moves between 
lattice sites, each atom in this molecule has to move one by one. This tunneling process 
naturally leads to the enhancement of molecular band mass in the strong coupling regime 
as M cx -EWci [isj] , where -Ebind is the molecular binding energy. Since -Ebind diverges in 
the strong coupling limit, the molecular mass M also diverges, leading to the vanishing 

. We note that the molecular mass 



2 



equals twice the atomic mass in the strong coupling regime of a uniform Fermi superfluid 



with no optical lattice, leading to the finite value of T c = 0.218Tp 
Fermi temperature). 



6, 



18 



20] (where Tp is the 



The Hubbard model is usually expected to be valid for the optical lattice system when the 
lattice potential is strong. However, Orso and co-workers [g, 9] recently studied a bound state 
problem in a realistic cosine-shape optical lattice potential, and showed that the molecular 
band mass M actually does not diverge but remains finite even in the strong coupling 
limit. Their result indicates that the Hubbard model is not valid at least for the strong 
coupling limit of optical lattice system. Thus, it is an interesting problem whether or not 
the molecular tunneling accompanied by dissociation predicted in the Hubbard model is 
really realized in a superfluid Fermi gas loaded on an optical lattice. This is also related to 
the problem about the validity of the Hubbard model for superfluid Fermi gases loaded on 
optical lattices. 

In this paper, we investigate a bound pair of Fermi atoms in an optical lattice. Including 
a cosine-shape periodic potential, we calculate molecular wavefunction. We show how the 
spatial structure of the bound pair changes during the tunneling through the lattice potential. 
We also compare molecular kinetic properties in the optical lattice potential with results in 
an attractive Fermi Hubbard model, in order to examine the validity of the Hubbard model 
for superfluid Fermi gases in optical lattices. 

This paper is organized as follows. In Sec. II, we explain a model optical lattice system, 
as well as how to calculate molecular wavefunction. Here, we also construct a single-band 
Hubbard model for the optical lattice system. In Sec. Ill, we calculate molecular excita- 
tions. They are compared with results obtained in the Hubbard model. In Sec. IV, we 
study the spatial structure of molecular wavefunction from the weak coupling regime to the 
strong coupling regime. We examine whether or not the virtual dissociation predicted in the 
Hubbard model occurs in the optical lattice system. In Sec. V, we consider the superfluid 
phase transition temperature T c in the BEC limit, where the Hubbard model gives T c = 0. 
Throughout this paper, we take h = ks = 1. We also set the system volume unity. 
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II. BOUND STATE IN MODEL ONE-DIMENSIONAL OPTICAL LATTICE 



We consider two attractively interacting Fermi atoms in a three dimensional system, in the 
presence of a one-dimensional optical lattice in the x-direction. These atoms are assumed to 
be in different hyperfme states, described by pseudospin a =|, j. The Hamiltonian is given 
byflfl 

H = # (r0 + H (t 2 ) - U5( ri - r 2 ), (2.1) 
where —U is the s-wave pairing interaction. The one-particle Hamiltonian Hq has the form 

t-t , n V 2 E r s / 2tcx\ , n 

where m is the mass of a Fermi atom. The last term in Eq. (12.21) describes the optical 
lattice in the x-direction, the height of which is measured in terms of the atomic recoil 
energy E r = 7r 2 /2md 2 . The lattice constant d is related to the wavelength A of laser light 
as d = A/20. In this paper, we ignore effects of a trap potential, for simplicity. 

Because of the contact pairing interaction in Eq. (12. ip . only the singlet pairing is allowed 
as the spin state of a bound pair. For the spatial part of the molecular wavefunction, noting 
that Eq. (12.11) is periodic in terms of the center of mass coordinate R = [x\ + x<i)j1 with 
the period d, one may take 

* q (rx,r 2 )= £ <C n2 (^pW r i)n+^(r 2 ). (2.3) 

p, 711,712 

Here, <pp(r) is an eigenfunction of the one-particle Hamiltonian Ho(r) in Eq. (12. 2p . with the 
energy e" where n is a band index. Since the system is uniform in the y- and ^-direction, 
the atomic energy has the form, = + (p 2 +p 2 z )/2m. Using the Bloch's theorem, we can 
write the eigenfunction in the form 0p(r) = e ip ' v Up x (x), where u" (x) is a periodic function 
satisfying Up x (x + d) = Up x (x). Because of the required antisymmetric property of fermion 
wavefunction, the spatial part \l/ q (r 1 ,r 2 ) must be symmetric with respect to the exchange 
of ri and r 2 . This is satisfied by imposing the condition <7"p m (q) = <?p 1,n2 (q) in Eq. (12. 3p . 

We note that (7p 1,n2 (q) with rii ^ n 2 describes interband coupling due to the spatial 
inhomogeneity by the lattice potential. In the extended zone scheme, this means that pairs 
of two atomic states with p and — p + G contribute to the molecular state when q = 0, 
where G is the reciprocal lattice vector. This is different from the case of a uniform gas, 
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where the molecular wavefunction with q = only involves pairs of atomic states with p 
and -p as ^ q ^ orm (r 1; r 2 ) = £ p g P e ip ri e- p r2 . 

Substituting Eq. (12.31) into the Schrodinger equation H^^ = E q ^/ q , we obtain 

n l5 n 2 / x = _ 

yp _m , n 2 771 

£ p+q/2 + £ -p+q/2 Q 

x E / ^ q/2 (rrn +q/2 (^ (2-4) 

k, 713, 714 

Since we are interested in the pair tunneling through the lattice potential, we take q = 
(q, 0, 0) in this paper. In addition, to examine the spatial structure of the molecular wave- 
function \l/ q (ri,r 2 ) in the ^-direction, we set y\ — yi = y and z\ = z 2 = z in Eq. (12.31) . 
The resulting molecular wavefunction ^(xx,^) = ^(^1, y, z\ x%, y, z) does not depend on 
y and z. Introducing the relative coordinate r = x\ — x 2 and the center of mass coordinate 
R — (x\ + X2)/2, one can rewrite ^^(xi, x 2 ) in the form 

^ q (r,i?) = ^ q (i? + r/2, R - r/2) 

= t MR ^ / -.-( q ) M - +g/2 (i ? + r /2)^ +9/2 ( J R-r/2). (2.5) 

The coefficient /^ 1,n2 (q) = S P!/ ,p z y'p 1 '™ 2 ^) obeys the equation 

/;, 1,n2 (q) = E - 



f "l _L_ F n 2 _ f] 

Py,Pz fc p+q/ 2 -P+q/2 ^ 



x E 5 / ^^/a(s)VW(s)^ 



fcai,>l3,™4 

In the cosine-shape periodic potential, (x) may be written as 



<(^) = E^"(0e^, (2.7) 



where C™ x (/) is determined by the equation 



Px , ^ 



+ — -e 



2m 2 



<£(0 - + 1) + - !)] = 0- (2-8) 



Substituting Eq. ( 12. 7ft into Eq. ( 12. 6ft . we execute the integration over a;. Then we obtain 

C n2 (q) = — -S — ztet E «(q;p«A)iE' n4 (q), ( 2 - 9 ) 

p v ,Pz £ p+q/2 £ -p+q/2 ^ k x ,n 3 ,n 4 

where 

«q;p*A) = E ^ 1 +9/2 (^)c'-L +? /2(^)^ 3 +? /2(^)^i +?/2 (i4). (2.10) 



The contact interaction in Eq. (12.11) brings about the ultraviolet divergence in Eq. (12.91) . 
Thus, we introduce the cutoff A c in the summations over (k x , n^, n^), as well as the cutoff Aj_ 



in the summation over pj_ = ^Jp 2 + p 2 in Eq. ( \2M . As usual, we eliminate effects of these 
momentum cutoffs by regularizing Eq. (12.91) . which is achieved by introducing the two-body 
scattering length a s given by 

Ana, U 



m 1 — all 

Here, a has the form 



(2.11) 



mA ±[ A c A 2 C + A 2 ± _ x A c] 

a= V V = —r^r -7— In . 9 + 2 tan — . (2.12) 

When we take A c /Aj_ ^> 1, Eq. (12.111) reduces to the familiar expression, 

Aira, U U 



(2.13) 



i-^E^ 1-17 2^" 
Executing the summations over p y and p z in Eq. (12.91) . one finds 

/;r 2 (q)=r^( P .) £ <;^(q;Px,^)^ 3 ' n4 (q), (2.14) 



where 

l. a2 

l px+q/2 C -px+g/2 ' 



l q VPx)- t— in 7 ni „ 2 , „ • 

We numerically solve the eigenvalue equation (12.141) to determine the molecular excitation 
spectrum E q , as well as /^'"^(q). The molecular wavefunction \l/ q (r, i?) is calculated from 
Eq. (12.51) . Since the current experiments on superfluid lattice Fermi gas are using a weak 
optical lattice potential[l[ Q], we take s = 3. For the momentum cutoffs A c and Aj_, it 
is difficult to take very large values because of computational problem. In this paper, we 
choose the value of A c so as to be able to include energy bands up to n = 14. For the 
cutoff Aj_, we set Aj_ = A c /6. Although numerical results on i? q still weakly depend on 
A c and Aj_, we expect that the essence of our results would be unaltered even when larger 
values of A c and Aj_ are used. We also find from numerical results that the truncation 
of the summations over (p x ,n 1 ,n 2 ) at A c in Eq. (12.51) affects the spatial structure of the 
molecular wavefunction around r = 0. For this problem, using the fact that one-particle 
wavefunction 0p(r) reduces to the plane wave in the high energy limit, we take into account 




FIG. 1: Calculated interaction Uu in the effective Hubbard model in Eq. (|2.16p . as a function of 
the pairing interaction measured in terms of the inverse scattering length a s . The inset shows the 
molecular binding energy -Ebind obtained from Eq. (|2. 14[) . 

the contribution from higher momentum region than A c by approximating 0™(r) to the plane 
wave in calculating \I/ q (r, R). Although this prescription cannot completely eliminate cutoff 
effects on \l/ q (r ~ 0, R), we can still study interesting molecular tunneling properties, using 
the spatial structure of \I / q (r, R). 

Besides the pair wavefunction \l/ q (r, i?), we also consider a single-band Fermi Hubbard 
model for the present periodic potential model in Eq. (12. ip . The single-band Hubbard model 
in momentum space is given by 



U 



u — t H V* t t 

n ~ 2^ £ P ( W ~~ "jy" 2^ C p+q/2T C -p+q/2| C -p'+q/2| C p'+q/2t 



p,cr 



(2.16) 



p,p',q 



where N is the total number of lattice sites in the x-direction. c' pa . is the creation operator of 



a Fermi atom with pseudospin a =],[. In Eq. (12.161) . we take the band dispersion e p so as 
to be eo.ua. tothe .owest energy band ea.eu.ated in the periodic potentia. system given 
by Eq. (12.21) [22]. The attractive interaction —Uh is taken so that Eq. fl 2 . 1 6 [) can reproduce 
the molecular binding energy -Ebmd = l-^q=o| obtained in the original periodic potential 
model given by Eq. (12.11) . Setting the molecular state as \^u) = J2 P ^p c P + q /2t c -p+q/2||0)' 
we obtain the equation for the energy of a bound state as 

U„^ 1 



i = ^j2 

N p £ P+q/2 + £ -p+q/2 



q 



Thus, Uh is given by 



Ui 1 = -Y 

N 



2e v + E } 



[2.17) 



(2.18) 



bind 
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FIG. 2: Molecular excitation spectrum AEq = £" q — Eq = Q. We take q = (g, 0, 0) and s = 3. 
AEq = — E^ =0 is the result obtained in the Hubbard model. q 2 /4m* is the molecular kinetic 
energy, when we assume a uniform system except that the atomic bare mass m is replaced by the 



band mass of the lowest atomic band, given by m* = (d 2 e p /dp 



Figure [T] shows the calculated Hubbard interaction Uh- In obtaining this result, we have 
used the binding energy -Ebind obtained from Eq. (12.141) . which is shown in the inset of Fig{TJ 
We briefly note that, in a periodic potential, a two-body bound state is possible even for 
negative scattering length 7, {J. (See the inset in FigJTJ) In a uniform system, a two-body 



bound state is possible only when a s 1 > 0. 
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III. EXCITATIONS AND BAND MASS OF A BOUND PAIR OF FERMI ATOMS 



Figure [2] shows the molecular excitation spectrum AE^ = — Eq=o (q = (g, 0, 0)). 
When the pairing interaction is not strong (panels (a) and (b)), we find that the Hubbard 
model in Eq. (12.161) well reproduces the excitation spectrum AE^ obtained in the periodic 
potential system. On the other hand, when (a s /c?) _1 = 2 (panel(c)), the Hubbard model 
underestimates AE^. 

Figure [3] shows the molecular (band) mass M in the x-direction, defined by 

M^r^r 1 . (3.i) 

^ dq 2 x /q-*> v ' 

In the weak-coupling regime {{a s /d)~ l ~ —1), the magnitude of the molecular mass M 
is close to twice the atomic band mass, given by m* = (d 2 e p /dpl)~\, . This can be also 
seen in Figj2](a), where the excitation spectrum AE q is well approximated to g 2 /4m* in the 
small momentum region. When (a s /c/) _1 ~ —1, since the molecule is weakly binding, the 
molecular motion is dominantly determined by the sum of two atomic band motions, which 
leads to M ~ 2m*. 

As one increases the strength of the pairing interaction, Figj3] shows that the molecular 
band mass M becomes heavier than 2m*. (See also Figsf5](b) and 2(c).) This mass enhance- 
ment can be described by the Hubbard model (M# in Figj3j) when (a s /c/) _1 < 1. However, 
when (a s /d)~ l > 1, the Hubbard model overestimates the molecular mass. In the strong 
coupling limit, while M approaches a constant value 8|, M# diverges 181 ] . 

The difference between M and Mjj in the strong coupling regime originates from different 
tunneling mechanisms between the original optical lattice system and the effective Hubbard 
model. As mentioned in the introduction, molecular motion in the Hubbard model is ac- 
companied by virtual dissociation in the strong coupling regime 18j . To see this in a simple 
manner, we consider the model shown in FiglU where a tightly bound molecule with the 
binding energy -Ebind moves from the i-th site to the (i + l)-th site. Noting that the creation 
of the intermediate state in Figfjjb) costs ffhipri , we obtain the nearest-neighbor molecular 
hopping matrix element as — tju = — 2t 2 IE^a [18|. where —t is the nearest-neighbor atomic 
transfer matrix element. When we only retain the tunneling process shown in FigJUby as- 
suming a small t, we obtain the molecular band = —2tMCOs(q x d), giving the molecular 
band mass M = E^nd/ {2td) 2 . Because -Ebind — ► oo in the strong-coupling limit, M diverges. 
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FIG. 3: Molecular band mass M = (d 2 Eq/dq 2 ,)~^ >Q , as a function of the inverse scattering length 
a s . Mh shows the result obtained in the Hubbard model, given by M# = (d 2 E^ / dq 2 )~^ , where 
Eq is determined by Eq. (|2.17p . M mo \ is the molecular mass obtained from Eq. () 3 . 3 j) . 





i+i 



FIG. 4: Simple model of pair tunneling in the strong coupling regime of the tight-binding model 
with the nearest-neighbor atomic hopping —t. Solid and open circles represent lattice sites and 
atoms, respectively. When a molecule moves from the i-th site to the (i + l)-th site, each atom has 
to move one by one. As a result, the modulate dissociates into two atoms in the intermediate state 
(panel (b)), so that the energy in the intermediate state is higher than that in the initial state by 
the binding energy i^bind- This tunneling process leads to the molecular hopping matrix element 
— tjtf = 2t 2 /E hm d- 

Although this is a simple evaluation, the enhancement of Mg in the strong coupling regime 
shown in Figj3] is found to be directly related to the molecular tunneling accompanied by 
virtual dissociation. 
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On the other hand, Orso and co-workers 



showed that the molecular tunneling mecha- 



nism in the strong coupling regime of cosine-shape periodic potential system is quite differ- 
ent. In this regime, since the molecular size is much smaller than the lattice spacing d, the 
r-dependence of the molecular wavefunction x I / q (r, R) is close to that in a uniform system, 
given by 

*(f) = -^=-e~ f ' a % (3.2) 



where r = ^Jx 2 + y 2 + z 2 . When we extract terms depending on the center of mass coordi- 
nate R from Eq. (12 .11) . we obtain 8] 

1 d 2 nr 2ixR 

Hn = — — — — E S S COS — — COS — 

H AmdR 2 d d 

1 d 2 2dE s s^ _ l7 ra s 2nR 

tan — — cos — - — . (3.3) 



4m dR 2 7ca s 2d d 

In obtaining the last expression, we have replaced the factor cos(7rr/d) by the expectation 
value (^(f)! cos(7rr/d)|^ r (f)) = (2d/na s ) tan _1 (7ra s /2(i). In the strong-coupling limit (a^ 1 — > 
+oo), Eq. (13.31) reduces to 

Equation (13.41) shows that a molecule feels the periodic potential with the finite height 
2E s s in the strong coupling limit js]. Namely, the molecular band mass M remains finite, 
in contrast to the case of Hubbard model. As shown in FigOl the molecular band mass 
M mo \ calculated from Eq. (13.31) explains the behavior of M in the strong coupling regime. 
This means that the bound pair moves without dissociation in the strong coupling regime of 
optical lattice system. 

Although the molecular tunneling accompanied by the dissociation does not occur in the 
strong coupling regime of the optical lattice system, we can still expect the possibility of this 
interesting tunneling phenomenon somewhere in the BCS-BEC crossover region, especially 
in the intermediate coupling regime ((a s /<i) _1 ~ 0), because the enhancement of M in this 
regime is in good agreement with Mh, as shown in Figj3j In the next section, we explore 
this possibility, based on the analysis of molecular wavefunction. 
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FIG. 5: Calculated molecular wavefunction ^ / q= o(r, i?) as a function of the relative coordinate 
r = x\ — x%. In this figure, we show |r x I , q = o(r - , R)\ to eliminate the r _1 behavior, which also 
appears in the molecular wavefunction in a uniform system. The coefficient C is chosen as C _1 = 



X^oo dr\r^fq(r, R)\ 2 . The dashed line in panel (a) shows ^ q= o(r) in Eq. (|4.ip . Sharp peaks and 
dips seen around r = may be artifacts, originating from the finite cutoff A c used in numerical 
calculations. (See the text.) 

IV. MOLECULAR WAVEFUNCTION AND VIRTUAL DISSOCIATION IN OP- 
TICAL LATTICE 

Figure [5] shows the calculated molecular wavefunction ^ q=0 (r, R), as a function of the 
relative coordinate r = x± — Xi. In this figure, the upper and lower panels show the cases 
when the molecular center of mass position is at the bottom of the periodic potential R = 
and at the top of the potential R/d = 0.5, respectively ((A) and (C) in Figj6j respectively). 
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V(x) 



FIG. 6: Center of mass position of the molecular wavefunction ^ q (r, i?) in Figj5j V(x) = 
-(E s s/2)[cos(2irx/d) - 1] is the periodic potential. A: R/d = 0. B: R/d = 0.25. C: R/d = 0.5. 
Results for D (R/d = 0.75) and E (R/d = 1) are the same as those at R/d = 0.25 and R/d = 0, 
respectively. 
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FIG. 7: Spatial variation of the molecular wavefunction ^(r, R) in the case of finite q = (q, 0, 0). 

The middle panels show the case of R/d = 0.25 (and also R/d = 0.75) ((B) and (D) in 
Figj6]). Namely, Fig. [5] shows how the spatial structure of the molecule varies when the 
molecule moves from (A) to (E) in Figj6j 

When the pairing interaction is weak (Figsj5](a)-(c)), the wavefunction spreads out. Ex- 
cept for panel (b), one finds oscillating structures, originating from the presence of optical 
lattice potential. For a given binding energy -Ebind an d atomic band mass m*, when we 
calculate the molecular wavefunction ignoring other lattice effects, we obtain 



^q=o(r) 



n 

_ p -Vm,*E hind \r\ 



(4.1) 



where C is a normalization constant 



2l| . Apart from the oscillating structure, the overall 



spatial variation of ^ q=0 (r) can be described by Eq. (14. ip . as shown in Figj5|a). This means 
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that the molecular size in this regime is dominated by the binding energy E^md an d atomic 
band mass m* appearing in Eq. (14.11) . 

Comparing FigJS^a) with Fig 0(c), we find that the peak positions in ^ q=0 (r) are different 
between the two. In panel (a), in addition to the central peak at r = 0, satellite peaks can 
be seen at r/d = ±2, ±4, • • •. In panel (c), the satellite peaks appear at r/d — ±1, ±3, • • •. 
In the former case, to satisfy R = avoiding the potential energy, one should put two atoms 
at (x\, X2) = (0, 0), (±d, =pi), (±2d, =F2<i), • • •. In the relative coordinate r = x\ — 22, these 
configurations immediately explain the peak positions in Figj5]^a). In the same way, the 
configurations which satisfy R/d = 0.5 and avoid the potential energy loss are (xi,X2) = 
(0, d), (d, 0), (—d, 2d), ■ ■ ■. These again explain the peak positions in FigJ5](c). The case of 
R/d = 0.25 is considered to involve both configurations, so that the oscillating structure is 
cancelled out to disappear, as shown in Figj5](b). 

As one increases the strength of the pairing interaction, the molecular wavefunction 
shrinks, reflecting the increase of the binding energy -Ebind- In the intermediate coupling 
regime shown in FigsJ5(d)-(f), while no satellite peak can be seen in panels (d) and (e), one 
finds two satellite peaks at r = ±d in panel (f). This means that the bound pair partially 
dissociates into two atoms at x = and x = d, when the center of mass position is at the 
top of the lattice potential R/d = 0.5. As shown in FigJTl these satellite peaks at r = ±d 
also exist when the molecule has a finite momentum q in the ^-direction. Thus, we find 
that the pair tunneling accompanied by dissociation predicted in the Hubbard model really 
occurs in the intermediate coupling regime of optical lattice system. 

However, in the strong coupling regime shown in FigsJ3(g)-(i), the satellite peaks are 
absent even when the molecular center of mass position is at R/d = 0.5 (panel (i)). Namely, 
the molecule tunnels through the lattice potential without dissociation in this regime, which 
is consistent with the discussion in the previous section. 



V. SUPERFLUID PHASE TRANSITION TEMPERATURE IN THE STRONG- 
COUPLING LIMIT OF LATTICE FERMI GAS 



Since the pair mass M in the optical lattice system actually remains finite in the strong 
coupling limit, we can expect a finite superfluid phase transition temperature T c even in the 



BEC limit, where the Hubbard model predicts the vanishing T c 
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19) . In this section, 
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FIG. 8: Dispersion of molecular excitation spectrum AE^ = E^ — E q= o in the (^-direction in the 
BEC limit (solid line). We take s = 3 and q = (q x ,0,0). This result is obtained from the model 
given by the Hamiltonian in Eq. (|5.3p . Solid circles show the result in the tight-binding model 
given by the first term in Eq. (|5.4p . where tM shown in the inset is used, tu is determined so that 
the band width Yltu of the Bose Hubbard model in Eq. (|5.4p can be equal to that of the lowest 
energy band obtained from the model Hamiltonian in Eq. (|5.3p . 

we evaluate T c in the BEC limit, including the finite value of M. 

For this purpose, we consider the BEC limit of a two-component Fermi gas in a three- 
dimensional cubic optical lattice. The Hamiltonian is given by 

H = J2H (r J )-Uj2S(r l ~r 1 ), (5.1) 

where the one-particle Hamiltonian H (r) has the form 

Jfo(r) = -— + — (3 - cos — - cos — - cos — ) . (5.2) 

In the BEC limit, Cooper pairs have been already formed above T c and the pair size is much 
smaller than the lattice constant d. In this case, one can treat the Cooper pairs as point 
bosons. Thus, instead of Eq. (15. ip . one may consider the simpler Hamiltonian, 

H = H Rx +H Ry +H Rz , (5.3) 

where H R is given by Eq. (13.41) . In addition, when the lattice potential is strong, the 
motion of bosons can be described by the tight-binding model with the nearest-neighbor 
boson hopping — tju- Indeed, as shown in FigJHl the molecular excitation spectrum obtained 
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from the model Hamiltonian in Eq. (15.31) is well approximated to the tight-binding dispersion 
= —2tM J2j= x ,y,z[ cos Qjd—1] when s > 3. We also note that, when the band gap between 
£p =1 and £p =2 in the original fermion system is very large in a strong lattice potential, one 
can ignore multiple occupation of bosons. These situations can be conveniently described 
by the Bose Hubbard model with infinitely large on-site repulsive interaction Um — ► +00, 

H M = -t„ £ [b% + h.c] + nf{nf - 1), (5.4) 

where b\ is the creation operator of a (molecular) boson, and nf 4 = b\bi. The first term 
describes the boson hopping between nearest-neighbor sites, where the summation is taken 
over the nearest-neighbor pairs. 

When the original fermion system is in the half-filling, the corresponding boson density 
equals % = 0.5 per lattice site. In this case, when we describe the occupied site and vacant 
site by pseudospin j and j, respectively, the Bose Hubbard model with Um — * +00 in Eq. 
( 15. 4p can be mapped onto the three-dimensional XY-model, 

h xy = t M j2( s i s -+ s i s -) 

(m) 

= 2t M J2(SiSi + SiSi). (5.5) 

Here, S± = S l x ± iS l y1 and S l x and 5j are S = 1/2 spin operators. Evaluating the phase 
transition temperature T c of this spin system within the simple mean-field approximation, 
one finds 

T c = 3t M . (5.6) 

We note that this result can be also obtained from Eq. ( 15.41) without mapping onto the spin 
model. We explain the outline of this alternative derivation in the Appendix. 

Figure [9] shows the calculated T c in the strong coupling BEC limit of a superfluid Fermi 
gas loaded on the cubic optical lattice (half-filling case). Comparing this result with the 
maximum T r ~ O.OAep obtained in the intermediate coupling regime of the Fermi Hubbard 



model 
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161 ]. we find that, although T c in the BEC limit remains finite due to the finite 



magnitude of M, it is still low when the lattice potential is strong (s ^> 1). 



16 



0.05 

0.04 
£j 0.03 
^ 0.02 

0.01 


4 8 12 16 20 
S 

FIG. 9: Superfluid phase transition temperature T c in the BEC limit of the superfluid Fermi gas 
loaded on the three-dimensional optical lattice, as a function of the potential height measured in 
terms of the atomic recoil energy E r . This figure shows the case of half- filling, e-p is the Fermi 
energy. 

VI. SUMMARY 

In this paper, we have investigated tunneling properties of a bound pair of Fermi atoms 
in an optical lattice. Including a realistic one-dimensional cosine-shape periodic potential, 
we have calculated the molecular wavefunction, binding energy, excitation spectrum, and 
band mass. We have also discussed validity of the Hubbard model for superfluid Fermi gases 
in optical lattices. 

In the strong coupling regime of the Hubbard model, the molecular tunneling through 
the lattice potential is accompanied by virtual dissociation into two atoms. This tunneling 
phenomenon does not actually occur in the strong coupling regime of real optical lattice 
system, where the bound pair moves in the lattice potential without dissociation. However, 
in the intermediate coupling regime, spatial structure of the molecular wavefunction indicates 
that the molecule dissociates into two atoms when the center of mass position is located at 
the top of the lattice potential. Our results show that the tunneling mechanism accompanied 
by virtual dissociation is realized in the optical lattice system, not in the strong coupling 
regime, but in the intermediate coupling regime. 

We have also examined the superfluid phase transition temperature T c in the strong 
coupling BEC limit, where the simple Fermi Hubbard model is no longer valid. Including 
the correct molecular tunneling process in this limit, we showed that T c is finite but is still 
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low compared with the maximum T c obtained in the intermediate coupling regime of the 
Hubbard model. 

So far, the BCS-BEC crossover in a lattice Fermi gas has been mainly examined based 
on the Hubbard model consisting of the atomic hopping and on-site pairing interaction. 
Since this model gives vanishing T c in the strong coupling limit, the observation of a finite 
and constant T c in the BEC regime of a lattice Fermi gas would be an evidence of the 
pair tunneling without dissociation. Although T c in this regime is expected to be low, the 
observation of the finite T c in the BEC regime is an important challenge to clarify the validity 
of the Hubbard model in considering the BCS-BEC crossover regime of lattice Fermi gases. 
We also note that the molecular tunneling with virtual dissociation enhances the molecular 
mass. This naturally leads to the suppression of T c in the intermediate coupling regime. 
Thus, the observation of the decrease of T c in the intermediate coupling regime would be an 
indirect evidence of the virtual dissociation of the bound pair during the tunneling through 
the lattice potential. Since the superfluid Fermi gas in an optical lattice is an important 
many-body system in both cold atom physics and condensed matter physics, we expect that 
our results would be useful for the study of basic properties of this interesting system. 
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APPENDIX A: ALTERNATIVE DERIVATION OF EQ. SSM 

In the superfluid phase, we take 6, = $ + 5bi in Eq. (15.41) . where $ is the BEC order 
parameter. In the mean field approximation, ignoring the fluctuation term having the form 
5b\5bj, we find that Eq. (15.41) reduces to the sum of the on-site Hamiltonian as Hu — 
T,iH M {i), where 

H M (i) = zt M ® 2 - zt M *(bi + 6j) - mn f - ^fnf{nf - 1). (Al) 

Here, we have added the chemical potential term —[Im^ to Eq. (lAlj) . z = 6 is the number 
of nearest-neighbor sites in the cubic lattice, and $ is taken to be real. When Um — » +oo, 
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one may only consider the vacuum state |0) and the single occupied state |1) = b\\0) 
Diagonalizing the on-site Hamiltonian HM{i), we obtain the eigenenergies as 

1 r 



E± = zt M ® 2 - ^ n M ± sjtfu + 4(^m$) 2 • (A2) 



The free energy per lattice site is given by 

1 



2' 







2 



F = zt M <$> 2 - -n M -T\n 2cosh^ v / y u| / + 4(zt M $) 2 . (A3) 



The superfluid order parameter $ is determined so as to minimize the free energy in Eq. 
(IA4[) . which gives 

n 2 M + A(zt M $y p 



tanh^^| f + 4(zt M $) 2 . (A4) 
zt M * 

The equation for T c is obtained by setting $ = in Eq. flA4j) . When % = 1/2, the equation 
for the number density of bosons is given by 

~ = a , (A5) 

which gives \im = 0. Substituting this result into Eq. ( 1A4I) with $ — > 0, we obtain Eq. 
(EE}. 



[1] J. K. Chin, D. E. Miller, Y. Liu, C. Stan, W. Setiawan, C. Sanner, K. Xu and W. Ketterle, 
Nature 443, 961 (2006). 

[2] D. E. Miller, J. K. Chin, C. A. Stan, Y. Liu, W. Setiawan, C. Sanner, and W. Ketterle, Phys. 
Rev. Lett. 99, 070402 (2007). 

[3] For reviews, see, L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford Univer- 
sity Press, NY., 2003) Chap. 16; M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. 
Sen, and U. Sen, Adv. in Phys. 56, 243 (2007). 

[4] E. Timmermans, K. Furuya, P. W. Milonni and A. K. Kerman, Phys. Lett. A 285, 228 (2001). 

[5] M. Holland, S. J. J. M. F. Kokkelmans, M. L. Chiofalo and R. Walser, Phys. Rev. Lett. 87, 
120406 (2001). 

[6] Y. Ohashi and A. Griffin, Phys. Rev. Lett. 89, 130402 (2002). 

[7] P. O. Fedichev, M. J. Bijlsma, and P. Zoller, Phys. Rev. Lett. 92, 080401 (2004). 

[8] G. Orso, L. P. Pitaevskii, S. Stringari, and M. Wouters, Phys. Rev. Lett. 95, 060402 (2005). 

19 



[9] G. Orso and M. Wouters, Phys. Rev. A 73, 012707 (2006). 
[10] R. B. Diener and T. L. Ho, Phys. Rev. Lett. 96, 010402 (2006). 

[11] M. Machida, S. Yamada, Y. Ohashi, and H. Matsumoto, Phys. Rev. A 74, 053621 (2006). 
[12] H. Tamaki, Y. Ohashi, K. Miyake, Phys. Rev. A 77, 063616 (2008). 

[13] A. Koga, T. Higashiyama, K. Inaba, S. Suga, and N. Kawakami, J. Phys. Soc. Jpn. 77, 073602 
(2008). 

[14] G. Watanabe, G. Orso, F. Dalfovo, L. P. Pitaevskii, and S. Stringari, cond-mat/0806.4713. 
[15] J. Mentink and S. Kokkelmans, cond-mat/0806.4784. 

[16] For a review, see, R. Micnas, J. Ranninger and S. Robaszkiewicz, Rev. Mod. Phys. 62, 113 
(1990). 

[17] M. Kohl, H. Moritz, T. Stoferle, K. Giinter, and T. Esslinger, Phys. Rev. Lett. 94, 080403 
(2005). 

[18] P. Nozieres and S. Schmitt-Rink, J. Low. Temp. Phys. 59, 195 (1985). 

[19] C. Chien, Y.He, Q. Chen, and K. Levin, Phys. Rev. A 77, 011602 (2008). 

[20] M. Randeria, in Bose-Einstein Condensation, edited by A. Griffin, D. W. Snoke and S. 

Stringari (Cambridge University Press, NY., 1995), p. 355. 
[21] In obtaining Eq. (|4.1|) . we have taken A c — > oo. Equation (|4.ip is consistent with Eq. (|3.2|) 

when -Ebind = l/m*a%. 

[22] When we take s = 3 in Eq. (|2.2p . the energy band s^ 1 is found to still deviate from the well- 
known band structure — 2t[cos(p x d) — 1] in the tight-binding model with the nearest-neighbor 
hopping — t. Thus, setting e p = £p =1 in Eq. (|2.16|) implicitly means that long range hopping 
terms in the x-direction are taken into account. We note that e^ 1 is well approximated to 
—2t[cos(p x d) — 1] when s > 10. 



20 



